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A simple proof of quantum adiabatic theorem is provided. Quantum adiabatic ap- 

proximation is divided into two kinds. For Hamiltonian H (t/T), a relation between 

Q ' the size of the error caused by quantum adiabatic approximation and the parameter 

■ T is given. 
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^ ■ I. INTRODUCTION 



^ I A quantum system is described by its Hamiltonian. When the Hamiltonian is time- 

s' \ independent and its spectral decomposition is known, we can easily get the dynamical evo- 
lution operator (DEO). However it is usually impossible to get an analytic expression of 
the DEO when the Hamiltonian varies with time. In a pioneering paper. Born and Fock 
considered a kind of time-dependent Hamiltonian Ht that has a special form Ht_= H (s) 
where s = t/T and the spectrum of if (s) consists of purely discrete eigenvalues [ij. Their 
result and the extended results [2I, 0, 0, 11] on more general H (s) are all called quantum 
i/-^ ! adiabatic theorem (QAT). For a history of QAT we refer to Ref. p]. 
O \ Recently a debate arised on the validity of application of quantum adiabatic theorem 

(QAT) or quantum adiabatic approximation (QAA) 0, 0, H, 0, Eol, [TH, One thing can 



be sure is that the widely used simple condition for QAA is actually insufficient 
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O ■ Since QAA and the related Berry phase [l3j have a wide application in many fields jl4l], it 
is valuable to find new conditions for the approximation. The recent discussion on QAT or 
QAA is also stimulated by quantum adiabatic computation [ist; q uantum computers are 
believed to be more powerful than computers that we use today 

To eliminate ambiguity, the thing should be specified first is the rigorous definitions 
of QAT and QAA. We only concern time-dependent matrix Hamiltonians without energy 
degeneracy, which are also the topics of the recent debate. QAT can only be discussed for 
Hamiltonians that have special form Ht = H (s) where s = t/T. This can be seen from the 
original papers on QAT [l|, 0] ■ We think QAT states that the difference between the DEO 
and the adiabatic evolution operator (AEO) (defined below) of the system will approach 
zero in the limit T — > 00. However, we can talk about QAA for Hamiltonians that may not 
be written in the form Hf = H (s). We define QAA as an approximation that uses the AEO 
in place of the DEO in our calculations. Under this definition of QAA, we find that there 
are two kinds of QAA and conditions for them to be acceptable are different. Of course we 
can also discuss QAA for Hamiltonians that are written in the form Hf = H{s), one can 
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discuss the relation between the amphtude of the error caused by QAA and the value of the 



parameter T [17 



The structure of the paper is as follows. In Sec. [TTl we give an introduction to QAT 
and present a proof of the theorem. In Sec. Illlt we give a discussion on QAA. In Sec. HVl 
we give an example to demonstrate the results obtained in Sec. IIIII In Sec. |Vl we give a 
relation between the size of the error caused by QAA and the parameter T for Hamiltonian 
H {t/T). Finally the conclusion is given in Sec. IVIi 



II. QUANTUM ADIABATIC THEOREM 

We just consider non-degenerate two- level systems to demonstrate our basic idea through 
out the paper. A smooth time-dependent Hamiltonian Hf has two instantaneous eigenvalues 
and eigenstates, 

Ht \mt) = mt \mt) , m = 1,2. (1) 

Eq. ([1)) cannot determine the phases of the eigenstates \mt). We choose the phases such 
that \m] 

{mt\ j^\mt) = 0, m= 1,2. (2) 
The DEO is denoted by Ud [t) which satisfies the Schrodinger equation 

'^^^ = HtUd (t) , Ud (0) = /. (3) 

Except for the trivial cases where the projectors \mt) {mt\ are independent of time t, it is 
usually impossible to obtain an analytic expression for Ud (t). The AEO is defined as 

2 

Ua{t)=J2(^-TJo-^"^^'\nt){no\, (4) 

n=l 

which is a little different from the so-called adiabatic transformation [3|. 

Theorem 1 (QAT) When Ht has the special form Ht = H (s) where s = t/T G [0, 1], i.e., 
it varies from H (0) to H (1) using time T , for any initial state |\E' (0)) there is 

hm ||[/,(t)|vI/(0))-f/,(t)|vI/(0))||=0. (5) 

1 — >oo 

Remark. What we consider is the simplest kind of QAT and it has been proved long 
ago by considering operator evolutions [l|, 0, H, 0], and no one has doubt on it. From the 
theorem it is not hard to see that when the system starts from the ground state of H (0) 
it will evolve closely to the ground state of H {!) if T is big enough; this is the basis of 
quantum adiabatic computation [l5| . Now we give a proof by considering the state vector 
evolution, which we think is more intuitive, and the process of the proof will be used when 
we discuss QAA. 

Proof. We write the system state at time t as 

2 

Ud (t) 1^ (0)) = 5^c„ {t) e-Ho*"*"^*' \nt) . (6) 

n=l 
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Since Ud (t) (0)) satisfies the Schrodinger equation, we can obtain [18 



dt 

dc2 (t) 
dt 



dt 

-c, [t) (2,\^\l^eifo<^^'-'^'^^r 
dt 



(7) 
(8) 



Integrate both sides of Eq. ([7]) we get 

ci (t")-ci (0) 



C2(t) (1,1 I |2,)eHo*(i.-2,K'rfi 



C2(t) 



i{it-2t)/h 



d 



(9) 



In the above we have used the fact that there is no energy degeneracy, i.e., 1, — 2, 7^ 0. 
Integrate by part we get 



where 



ci it") - ci (0) = A it") + B it") + C it") 
Ait") = -C2 it) fteTJoih'~2„)dt' 



t=t" 

t=0 ' 



B it") 
C it") - 



c, it) |ieno(V-2,)^*'rft, 



dt 



ft 



dt 

(1*1 ||2*) 



dt, 



(10) 

(11) 
(12) 

(13) 

(14) 



iiit-2t)/h- 

Because H (s) depends only on the parameter s, instantaneous eigenstates and eigenvalues 
are also dependent only on s, i.e., we can write \mt) = \m is)), rrit = m (s). Now we have 



ft 



From Eqs. flllHlSI) we can get 



1 ^(1 is)\i\2is)) _ fis) 
T tilis)-2is)) T ' 



|A(m<-max|/(.) 

T se[o,i] 



(15) 



\B it")\ < - me.x 
T se[o,i] 



|C(ni<;^max 

I s6[0,l] 



pis] 
ds 

d 



(2 (3)1 ^|1 (.))/(.) 



(16) 

(17) 
(18) 
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In deriving Eq. (fT8|) from Eq. (fT3ll we have used Eq. ([8]). From Eq. (fTOl) . we know 
I A + 1 5 + |C {t")\ is an upper bound of |ci (t") - ci (0)|, and Eqs. ^MB indicate 
this upper bound will approach zero in the limit T — oo, so we have 



Similarly we can obtain 
Because 



we get 

which completes the proof. 



hm |ci(t)-ci (0)1 = 0. 

1 -^oo 



hm |c2(t)-C2 (0)1=0. 

I — ^■oo 



|[/,(t)|vi/(0))-t/,(t)|vI/(0))|| 



\ n=l 



hm \\U,{t)\^m-Ua{t)\<iI 

1 -^oo 



0, 



(19) 
(20) 



(21) 



(22) 



III. QUANTUM ADIABATIC APPROXIMATION 



The system state {"^ (t)) at time t and the initial state |\I/(0)) are connected by the 
relation (t)) = Ud (t) |\E' (0)). Suppose the operator B (t) represents a physical quantity 
at time t. Quantum mechanics tells us that when we measure the quantity the result will 
be random and the average value will be 

{t)\B (t) \^ (t)) = (VI/ (0)1 f/] (t) B (t) Ud (t) \^ (0)) . (23) 

Assume we are given the initial state |^ (0)) and the Hamiltonian Hf of the system. If we can- 
not figure out Ud (t) from H^, usually we wih not know the average value (^^ (t) \ B (t)]"^ (t)). 
From QAT we know that in some cases the difference between Ud [t) and Ua (t) will be very 
small, we may consider whether it is acceptable to use Ua (t) in place of Ud (t) in calculating 
the average in Eq. fl23l) . We define QAA as an approximation that uses the AEO Ua (t) 
in place of the DEO Ud {t) in our calculations. The definition of QAA applies to a general 
time-dependent Ht contrast to QAT. If we require the approximation to be acceptable for 
any physical quantity B (t) and any initial state I^E' (0)) in Eq. fl2^ . it means the difference 
between Ud (t) |^ (0)) and Ua (t) |^ (0)) should be smaU. 

The difference between Ud (t) I^E' (0)) and Ua (t) |^ (0)) is small means 

||f/d(t)|^(0))-f/„(t)|^(0))||<l, (24) 

which is equivalent to 

(vl/(0)|f/t(t)f/,(t)|vl/(0))^l. (25) 

Due to the linearity of quantum mechanics, for arbitrary initial state I^P (0)), condition fl25|) 
will be satisfied when 

(mo I Ul (t) Ud (t) |mo) ^ 1, m = 1, 2, (26) 
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which we regard as the condition for the first kind of QAA. This kind of QAA pays attention 
to the phase of 11^ (t) |mo). Condition (126!) ensures that the relative phase between the two 
instantaneous eigenstates in Ud (t) (0)) is almost the same as that in Ua (t) (0)). 

The probability of finding the system in an instantaneous energy eigenstate \nt) is of con- 
siderable interest, e.g., in coherent population transfer among quantum states of atoms and 
molecules Il9|. In this case, in Eq. ( |251) the operator B (t) = \nt) {nt\ and the probability is 



I {nt\ Ud (t) 1^ (0)) 1^ From the definition of the Ua (t) in Eq. we know | (rij | Ua (t) |^ (0)) |^ 
is a constant, it means the system will follow the instantaneous eigenstate if it starts from 
an instantaneous eigenstates. Therefore it is acceptable to use Ua (t) in place of Ud (t) in 
calculating the probability \{nt\ Ud (t) \mQ)f when 



Eq. ( I27I) is the same as 



\{mt\ Ud{t) |mo)| ^1, m = 1,2. 



|(mo|?7i(t)f/rf(t)|mo)| ^1, m = 1,2, 



(27) 



(28) 



which we regard as the condition for the second kind of QAA. This kind of QAA is considered 
in ji, 18|. When condition (l28l) is satisfied, the relative phase between the two instantaneous 
eigenstates in Ud [t] \^ (0)) may differ much from the relative phase in Ua {t) |^ (0)). 

The second kind of QAA is just a special case of the first kind and the conditions for 
them are different: condition fl26|) can lead to fl28|) while fl28|) may not lead to fl26|) . so 
there are cases where the second kind of QAA is acceptable while the first kind of QAA is 
inacceptable. Though the conditions for two kinds of QAA are given, usually it is not easy 
to directly check whether they are satisfied. In the following we will give a discussion on the 
conditions for them. 

From (12 ip we know that if 

|c„(t) -c„(0)| < 1, n = 1,2, 
the first kind of QAA will be acceptable. From Eq. (fTOl) we can get 

ici (n-ci m<\A{t")\ + \B{t")\ + \c{t")\, 



where 



\A{t")\ = 
\B{t"] 
\C{t")\ 



Mi*l||2*) 



1.-2. 



t max 

*G[0,t"] 



+ 



i=0 



t=t" 



d {^{U\i\2t) 



dt 



t max 

te[o,t"] 



lt-2t 



(29) 
(30) 

(31) 
(32) 
(33) 



Eq. (1501) gives an upper bound for |ci {t") — C\ (0)|. The first kind of QAA will be acceptable 
if |v4(t")|, \B{t")\ and \C are all very small. Eqs. ([32]) and ([33]) indicate that it 
may be inappropriate to use Ua (t) {"^ (0)) in place of Ud {t) |\I' (0)) when t is large. When 
instantaneous eigenstates and eigenvalues are given, it may be easy to check whether \A [t") | , 
\B (t")| and |C* (t")| are small at the given time t". For the example we will discuss in the 
next section, this method is quite good. 
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The second kind of QAA is discussed in Refs. Il8|. It is known that when 



d 



dti 



<t: 1, n ^ m, 



(34) 



the second kind of QAA will be acceptable [1, [l8|. We emphasize that the instantaneous 
eigenstates in (IMIl satisfy the phase condition ([2])_. We can derive Eq. flMl) as follows. 
It is not hard to know that the unitary operator U (t) = Ul (t) Ud (t) is generated by the 
Hamiltonian 

a 
di 



-thJ2 {ml ^ In*) e^/o("^*'-"*')'^*' |mo) ( 



no . 



(35) 



We can get an expansion of U (t) 

1 



U{t) = I + 



iH 



* ~ 1 
dtiHt, + 



thy 



dti / dt2Ht,Ht, + 



(36) 



Substitute Eq. (jlj) and Eq. ( l36l) into the expression Ud (t) = Ua (t) U (t), we can obtain an 
expansion of Udif). Similar discussions appear in Refs. |2o|, HH. The condition ( l28l) of the 
second kind of QAA is equivalent to 



\{mo\Ul{t)Ud{t)\no) \ = {mo\Uit)\no) 



<^ 1, m ^ n. 



(37) 



When we just substitute the first two terms of Eq. (136|) into (!37|) . we will obtain the condition 
(|34l) . Though ( !34l) as a condition for the second kind of QAA has not been strictly proved 
since we just use the first two terms in Eq. we believe it is sufficient. If (m^J ^ |nj^) 

and {rriti — nf ) are constants, condition fl51|) can be simplified into 



h{mt,\ ^ \nt,) 



^ 1, m^n, 



(38) 



while when (m^J ^ |nt^) and (m^/ — rij') are not constants, a modification of (138!) such as 



max 

tie[o,t] 



^ 1, m^n, 



(39) 



cannot replace (IMI) as a sufficient condition for the second of QAA [6|, [lO, ll8|, 122] , however it 
may be acceptable to regard fl39|) as a sufficient condition when (mi^ | ^ jritj) and (mt/ — nti) 
vary slowly. The counterexample in Ref. [g] demonstrates that a contradiction will appear 
at a special evolution time if we take (!39|) as a sufficient condition. We think the rapid 
changing of (m^J ^ In^^) causes the contradiction because when (m^J ^ varies fast, 
especially when its varying frequency resonant to the energy gap, condition fl5^ deviates 
much from (IMj) . Finally we want to emphasize that even when condition (IMj) is satisfied, it 
doesnot means it is acceptable to use Ua it) |no), condition to approximate Ud if) |no). 
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IV. AN EXAMPLE 



In this section we will give an example to demonstrate that there are cases where the 
second kind of QAA is acceptable while the first kind is inacceptable. 
Consider the Hamiltonian 

Ht = —hjjQ [ax cos 2ujt + ay sin 2ujt] , (40) 

where ojq and uj are both positive. Its eigenvalues are It = huo and 2f = —Jtwq. We choose 
the eigenstates 

\lt) = (e-*"*|0)-e*"*|l))/v^, (41) 
\2t) = (e-''"*|0) + e'"*|l))/y2, 

which satisfy the phase condition ([2]). It can be proved that in the basis {|lo) , |2o)} 

j*"^"* (cos cut - ^ sin cot) e^'^o* x ^ sin cut 

^-tujot X gjj^ ^-lujot /(.Qg tut + ^ sin Ljt) 

UJ \ UJ J 



ul it) U, it) 



(42) 



where uj = ^/uJ^~+^J^. 

The condition (!28|) for the second kind of QAA can be written as 



{cos iJtf + sin Lut^ ^1, (43) 

which can be simplified into 

-l^sin'^LDt^O. (44) 

From (jUj) it can be seen that when uj <^ ujq the second kind of QAA will be acceptable no 
matter how lage t is. This result can also be obtained from Eq. fl34l) . 
The condition fl26|) for the first kind of QAA can be written as 

e*"^"* (^cos Ljt - ^ sina;t^ ^ 1, (45) 

which is the same as ^ 

e"^^ + ze^'^ot ^^jsincjt ^ 1. (46) 



Notice that (145|) can lead to (1431) . it is not hard to see that (146)) is equivalent to (I44j) plus 
uH/ {lD + ujq) ^ 1. When u <ti ujq and uH/ujq ^ 1 the first kind of QAA will be acceptable. 
Though condition c<j -C cuq is enough for the second kind of QAA to be acceptable, it cannot 
certainly make the first kind of QAA acceptable when t is large and this result has been 
implied in Ref. [23]. 

In this example, the upper bound expression Eq. (I5U]) can be written as 

|ci(t)-ci(0)|<- + ^. (47) 

From (ITTl) we can also conclude that when u uq and uH/uq -C 1 the first kind of QAA 
will be acceptable, which coincides with the result obtained by figuring out Ul (t) Ud (t). 
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V. QUANTUM ADIABATIC APPROXIMATION FOR HAMILTONIAN H [t/T) 



The above discussion on QAA applies to a general time- dependent Ht- When Ht has the 
special form Ht = H (s) where s = t/T G [0,1], we can discuss the relation between the 
parameter T and the error caused by QAA. 



Theorem 2 Suppose Ht has the special form Ht = H (s) where s = t/T G [0, 1]. When 



I > — — max 

se[o,i] 



dH{s) 


7 

- + — 


dH{s) 


2 




ds 

o 


ds 


- + - 


ds^ 



\ 



\ 



J 



A{sy A{sy A(s) 
where A (s) = |1 (s) — 2 {s)\, the error caused by QAA cannot be bigger than 6, i.e., 



(48) 



\\u,{t)\^>m-Ua{t)\^m\\<5- (49) 

Proof. In Sec. HI we write Ud if) |* (0)) in Eq. ^ and we know from Eqs. ffTITOfTSl) that 

(50) 



\Cn {t)-Cn (0)1 



< — max 

T se[o,i] 



where 



2 1/(^)1 + 



f{s) 



df{s) 



ds 



+ 



h{Hs)\i\2is)) 

t {1 is) - 2 is)) 



(51) 



is defined in Eq. f lT5]) . Now we give an upper bound of I Cn (t) ~ Cn (0)| expressed in norm of 
the Hamiltonian H (s) and its derivatives. First we have [4, [l8 



{lis)\-\2is)) 
ds 



{lis)\^ 



ds 



2 is)) 



2(s)-l(s] 



so there are 



and 



1/(^)1 



Ml(3)|^|2(s)) 



il is) - 2 is))^ 



h 



< 



dH{s) 
ds 



,2 ' 



(2(^)|^|l(^))/(^) 



< 



dH{s) 



ds 



It can be proved that 



ds 



d 
ds 



ds 



AisY 

Hs)) 



(52) 



(53) 



(54) 



-^ilis)-2is)y 



(55) 



h 



< 



d^H{s 



ds^ 



h 



dH{s) 
ds 
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Substitute Eqs. §MB into Eq. ([50D, we get 



h 

< — max 

T sG[0,l] 



dH(s) 
ds 



d^H(s) 



ds'^ 



dH(s) 
ds 



A(s] 



A is] 



A is) 



Because 



\UAt)\^m-Ua{t)\^{0))\\ 



\ n=l 



when 



HI) is satisfied we have 

\\U,{t)\^{0))-Ua{t)\^{0))\ 



< max 

T sg[o,i] 



<5, 



dH{s) 



ds 



d^H(s) 



ds'^ 



dH{s) 
ds 



A(s] 



A{sy 



A is) 



(56) 



(57) 



(58) 



(59) 



which completes the proof. 

A resuh similar to theorem [2] appears in Ref . 17|] , which is derived from a very different 
method. Under the same error 6, the required evolution time T in our result is proportional 
to 1/(5 while in Ref. it is proportional to 1/5^. We think this is one of the most differences 
between the two results. ■ 



VI. CONCLUSION 

We give a discussion on the conditions of QAA. We present a proof of QAT, which we 
think is easier than that appears in the textbook [4]- We think there are two kinds of QAA, 
one cares the relative phase in the approximate system state while the other does not. For 
the kind of Hamiltonian H {t/T) we give a relation between the size of the error caused by 
QAA and the parameter T. 
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